Lifshitz Point in the Phase Diagram of Resonantly Interacting 6 Li- 40 K Mixtures 
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We consider a strongly interacting 6 Li- 40 K mixture, which is imbalanced both in the masses and 
the densities of the two fermionic species. At present, it is the experimentalist's favorite for reaching 
the superfluid regime. We construct an effective thermodynamic potential that leads to excellent 
agreement with Monte Carlo results for the normal state. We use it to determine the universal 
phase diagram of the mixture in the unitarity limit, where we find, in contrast to the mass-balanced 
case, the presence of a Lifshitz point. This point is characterized by the effective mass of the Cooper 
pairs becoming negative, which signals an instability towards a supersolid phase. 

PACS numbers: 03.75.-b, 67.40.-w, 39.25.+k 



Introduction. — Ultracold atomic Fermi mixtures are 
presently at the center of attention of both experimen- 
tal and theoretical physicists. Due to the amazing control 
over this system many fundamental discoveries have been 
made, while more are likely to follow soon. These dis- 
coveries started with reaching the so-called BEC regime 
where the balanced two-component Fermi mixture turns 
superfluid due to the Bose-Einstein condensation (BEC) 
of molecules [E 0] ■ By using a Feshbach resonance to vary 
the interaction strength between atoms in different spin 
states, the BEC-BCS crossover between a Bose-Einstein 
condensate of molecules and a BCS superfluid of Cooper 
pairs could be directly observed [3( . Since pairing is opti- 
mal for an equal amount of atoms in each spin state and 
is absent for the noninteracting fully polarized system, a 
phase transition occurs as a function of spin imbalance 
P, [1| . The phase diagram of the polarized mixture was 
for strong interactions found to be governed by a tri- 
critical point that resulted in the observation of phase 
separation [f| Q . The presence of gapless Sarma super- 
fluidity is also likely to be present in this system [J Q , 
but has not been unambiguously identified yet. 

Most recently, experiments have indicated that the 
physical consequences of yet another parameter can be 
explored, namely that of a mass imbalance between the 
two fermionic components. A very promising mixture in 
this respect consists of 6 Li and 40 K, which has a mass 
ratio of 6.7. Several accessible Feshbach resonances are 
identified in the mixture Q , while both species have also 
been simultaneously cooled into the degenerate regime 
[Io| . So far, experimental interest has focused on the 
BEC side of the Feshbach resonance, where molecules 
are formed. Although the mass imbalance itself does 
not lead to fundamental new physics here, this situation 
changes when the heteronuclear molecules are optically 



pumped to their ground state Then the molecules 
acquire a large electric dipole moment, which gives rise 
to anisotropic long-ranged interactions. In an optical lat- 
tice, this can lead to supersolid phases [12]. 

In this Letter we focus on a different regime, namely 
the so-called unitarity limit, where the s-wave scattering 



length of the interspecies interaction diverges. Here, the 
size of the Cooper pairs is comparable to the average in- 
terparticle distance and the pairing is a many-body effect. 
The mass imbalance has a profound effect on the pairing 
now, because it strongly alters the Fermi spheres. We 
show that for the sufficiently large mass ratio of the 6 Li- 
K mixture, the phase diagram not only encompasses 
all the exciting physics known from the mass-balanced 
case, but is even much richer. Similar to the solely spin- 
imbalanced case is the presence of phase separation [13| , 
which can occur due to the mismatch of the Fermi sur- 
faces. Also similar is that gapless Sarma superfluidity is 
unstable at zero temperature [HI], while there is a pre- 
dicted crossover to the Sarma phase at nonzero temper- 
atures 0]. However, the most exciting difference that 
we find is the presence of a Lifshitz point in the phase 
diagram. 

At a Lifshitz point the transition to the superfluid 
phase undergoes a dramatic change of character. Rather 
than preferring a homogeneous order parameter, the sys- 
tem now becomes an inhomogeneous superfluid. This 
exotic possibility was early investigated for the weakly in- 
teracting mass-balanced case by Larkin and Ovchinnikov 
(LO), who considered a superfluid with a single standing- 
wave order parameter [l6| . This is energetically more fa- 
vorable than the plane-wave case studied by Fulde and 
Ferrell (FF) 0. Since the LO phase results in periodic 
modulations of the particle densities, it is a supersolid 
[3]. The FF and LO phases have intrigued the physics 
community for many decades, but so far remained elusive 
in experiments with atomic Fermi mixtures. Typically, 
Lifshitz points are predicted at weak interactions where 
the critical temperatures are very low. However, in this 
Letter we show that the very special phase diagram of the 
6 Li- 40 K mixture contains both a Lifshitz and a tricritical 
point in the unitarity limit, as shown in Fig. 1. This 
is in sharp contrast to the mass-balanced case, where at 
unitarity a large body of theory only finds a tricritical 
point, in agreement with experiments jfj. 

In first instance, all these expectations follow from 
mean-field theory, which is useful for a qualitative de- 
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FIG. 1: Universal phase diagram for the homogeneous 6 Li- 40 K mixture in the unitarity limit as a function of temperature T 
and polarization P. The temperature is scaled with the reduced Fermi temperature &bTf = ef = h 2 (3n 2 n) 2 ^ 3 /2m, where m is 
twice the reduced mass and n is the total particle density. The result of mean- field theory is shown in panel a). For a majority 
of light 6 Li atoms there is a tricritical point (TCP), at which the normal state (N), the homogeneous superfluid state (SF), and 
the forbidden region (FR) meet each other. For a majority of heavy 40 K atoms there is a Lifshitz point (LP), where there is 
an instability towards supersolidity (SS). The size of the supersolid stability region is not calculated within our theory and the 
dashed lines are therefore only guides to the eye. The grey area sets the scale for panel b), where fluctuation effects are taken 
into account to calculate the location of the Lifshitz and the tricritical point more accurately. 



scription of the physics. However, it cannot be trusted 
quantitatively in the unitary regime, since it vastly over- 
estimates critical temperatures in that case. It is thus 
important to understand how much these are lowered 
by fluctuation effects. In the mass-balanced case two 
effects are dominant, namely the fermionic selfenergies 
and the screening of the interaction due to particle-hole 
fluctuations 14[. Taking these into account gives in an 
intuitive manner results that compare well with Monte 
Carlo calculations [IH and experiment We find that 
an extension of this procedure to the Li- K mixture 
leads to a reduction of the mean-field critical temper- 
atures by an experimentally relevant factor of 3. This 
makes it more difficult, but not impossible, to reach the 
superfluid regime. Moreover, the presence of the Lifshitz 
point in the phase diagram is unaffected by the fluctua- 
tions, allowing for an experimental study of supersolidity 
in the 6 Li- 40 K mixture. 

Phase diagram. — The critical properties of the su- 
perfluid transition in a fermionic mixture are determined 
by the effective Landau theory for the superfluid order 
parameter A(x). Although we consider no external po- 
tential, the order parameter may still vary in space due to 
a spontaneous breaking of translational symmetry. Close 
to the continuous superfluid transition, we expand the 
exact effective thermodynamic potential as 



0[A1 = 



J dx| 7 |VA| 



0. 



, (1) 



where the challenge is to express the expansion parame- 
ters in terms of the temperature T and the atomic chem- 
ical potentials fi±. The upper (lower) sign refers to the 



light 6 Li (heavy 40 K) atoms. A phase transition has oc- 
curred when the global minimum of ft is located at a 
nonzero order parameter (A(x)), which describes a con- 
densate of pairs. When j(T) is positive, the pairs have a 
positive effective mass and their center-of-mass state of 
lowest energy is at zero momentum. Then, we can con- 
sider a homogeneous pairing field A, for which a second- 
order transition occurs at a critical temperature T c de- 
termined by a(T c ) = 0. 

A continuous transition only occurs when the mini- 
mum at small values of (A) is global, which is not neces- 
sarily the case. The expansion of f2 may contain higher 
powers of |A| 2 that have negative coefficients, leading to 
a first-order transition with a jump in the order param- 
eter when ft [0] = f2[(A)]. Second-order behavior turns 
into first-order behavior when (3{T) becomes negative, so 
that the temperature T c3 at the tricritical point (TCP) 
is determined by a(T c3 ) = and (3(T c3 ) = 0. Another 
intriguing possibility is that not f3(T), but rather j(T) 
goes to zero. This leads to a Lifshitz point (LP), which 
is thus determined by a{Ti) = and 7(Tl) = 0. Since 
the effective mass of the Cooper pairs becomes negative 
below the Lifshitz point, it is energetically favorable for 
them to have kinetic energy and form a superfluid at 
nonzero momentum. This can be established in many 
ways, namely through a standing wave II 61 or a more 
complicated superposition of plane waves [lil. [20j | . Due 
to the variety of possibilities it is hard to predict which 
lattice structure is most favorable. However, the fact 
that they all emerge from the Lifhitz point facilitates the 
experimental search for supersolidity in the 6 Li- 40 K mix- 
ture. Moreover, a nonzero gap (A) gives rise to a sizeable 
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superfluid fraction, showing that the LO-like incommen- 
surate supersolid is very different from the commensurate 
supersolid studied in 4 Hc. 

We continue our discussion at the mean-field level, 
which is useful for further explaining the relevant physics. 
Starting point is the BCS thermodynamic potential den- 
sity for the 6 Li- 40 K mixture with masses m± [l3j], 



^bcs[A;/v] = 



rfk f /, N * /, N l A | 2 

— { e(k) - /* - foj(k) ^ J - L 



(2tt) 

k B T In ( 

cr=± 



1 



2e(k) 



where we introduced the average chemical potential 
/i = (/x + + /i-)/2, the difference ft = fu+ — 
and the reduced kinetic energy e(k) = ?i 2 k 2 /2m with 
m = 2m + rri-/(m + +m_). The Bogoliubov quasiparticle 
dispersions become ?kiv(k) = ftw(k) — ct[2/i — e+(k) + 
e_(k)]/2 with 7iu;(k) = ^/[ e (k) - m] 2 + |A| 2 and e CT (k) = 
7i k 2 /2m ff . We can now apply the exact critical condi- 
tions to ojbcs and obtain the mean-field phase diagram. 
Although the BCS potential neglects fluctuations in the 
order parameter, it is expected that fluctuation effects 
only result in quantitative corrections. This expectation 
stems from the strongly interacting experiments for the 
mass-balanced case, for which the mean-field diagram 
is topologically correct 0. The coefficients determin- 
ing the second-order phase transition and the tricritical 
point are readily calculated as a — c^bcs [0; A| 2 
and (i — <9 2 cjbcs[0; ^ a ]/d 2 \A\ 2 . The results are shown 
in the phase diagram of Fig. la. The polarization is de- 
fined as P = {n + — H- ) / (n + + n_ ) , while the particle 
densities are determined by n a = — <9o->bcs[(A); /i,j]/9/i CT . 
Therefore, the polarization is discontinuous simultane- 
ously with the order parameter, which gives rise to a 
forbidden region (FR) below the tricritical point. 

From Fig. la, we see that the mean-field phase di- 
agram also contains a Lifshitz point. It is calculated 
from the noninteracting Green's function for the Cooper 
pairs G/\(k, iL0 n ). In the normal state, this propagator 
is given by hG^ 1 (k, ioj n ) = 1/T 2B — hE(k, ioj n ), where 
rp2B _ 4 na fi 2 j m w jth a the diverging scattering length, 
while 5 is the so-called ladder diagram, given by 



KB(k, iu/ v 



dk' f 1 

(2tt) 3 \ 2e(k') 
1 — N + (k') — iV_(k • 



(3) 



k') 



ihujn - e+(k') - e^(k - k') + 2^ 



with N a (k) = i/[ e (^(k)-M.)/fc B T + ]_] the Fermi dig _ 
tributions. Note that the mean- field expression for 
a is equal to —hG^ 1 (0,0), while we have that 7 = 
-dTiG~^(O,0)/dk 2 = at the Lifshitz point. What 
precisely happens below the Lifshitz point is an intrigu- 
ing question for further research. In Fig. la, we have 
sketched a simple scenario, where there is a second-order 



transition from the normal to the supersolid phase, for 
which the condition is G^ 1 (kss,0) = with kss the 
wavevector of the supersolid. However, this transition 
can in principle be of first order, where the realized super- 
solid periodicity can also contain more than one wavevec- 
tor [19( | . Moreover, the transition from supersolidity to 
the homogeneous superfluid phase is also expected to be 
of first order. The calculation for the stability regions of 
all possible supersolid lattices is beyond the scope of this 
Letter, such that the dashed lines in Fig. la are merely 
guides to the eye. 

Strong interactions. — Having established that within 
mean-field theory a tricritical point dominates the phase 
diagram when there is an abundance of light 6 Li atoms, 
while a Lifshitz points is present when there is an abun- 
dance of heavy 40 K atoms, the question arises whether 
this interesting physics remains when we take fluctuation 
effects into account. Answering this question also leads 
to more quantitative predictions for future experiments. 
The normal equation of state at unitarity is known to be 
strongly affected by fluctuations [2l| . Since this equation 
of state exactly follows from the fermionic selfenergies, we 
first try to take these accurately into account. We achieve 
this through a simple construction that gives excellent 
agreement with Monte Carlo results, and for which we 
only need to know the sclfcncrgy KZi a of a single a atom 
in a sea of —a atoms. At zero temperature, it is given 
by TiYi a = — Ccr/i-o-. The coefficients can be calculated 
within the ladder approximation 22] giving c+ = 2.2 and 
c_ = 0.44, or with a renormalization-group approach 14 1 
giving c + = 2.0 and c_ = 0.34, while Monte Carlo calcu- 
lations lead to c + = 2.3 and c_ = 0.36 Noting that 
these approaches agree rather well, we use the Monte 




FIG. 2: (Color online) Equations of state for the zero- 
temperature normal state of the 6 Li- 40 K mixture (full line) 
and the mass-balanced mixture (dashed line) in the unitar- 
ity limit. The Monte Carlo results for the 6 Li- 40 K mixture 
by Gezerlis et al. [23| are shown by circles, while the Monte 
Carlo results for the mass-balanced mixture by Lobo et al. 
[2^ ] are shown by squares. The energy density e is scaled 
with the ideal gas result efg = 3(eF+n+ + £F_rz_)/5 with 
£Fct = ft 2 (67r 2 n CT ) 2,/3 /2?Ti (T . The deviations from one thus show 
the strong interaction effects. 
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Carlo results for the rest of the calculation, since they 
effectively incorporate all Feynman diagrams. Moreover, 
KT,_ (7 = because the majority atoms are unaffected by 
the single minority atom. The chemical potential of the 
minority particle is then given by \i a = h~E a . 

Next, we define renormalized chemical potentials as 
fi' a = (i a + Cff/Z-oVO"^ + M-cr) 0- Two important solu- 
tions to these equations are such that the chemical poten- 
tial of one species (—a) is not renormalized, \J_ a = 
while the renormalized chemical potential of the other 
species is zero, i.e., \i! a = + c a n- a = 0. These two 
solutions correspond precisely to the two extremely im- 
balanced limits. By using the renormalized chemical po- 
tentials in the mean-field thermodynamic potential, we 
can calculate the full normal equation of state by using 
the equation for the densities n a — — <9(jJbcs[0; d[i a . 
The comparison with the Monte Carlo equation of state 
is shown in Fig. 2 and the agreement is excellent. We also 
show the comparison for the mass-balanced case, where 
both ladder and Monte Carlo calculations give c± = 0.61. 
The agreement shows that our construction captures the 
full thermodynamics of the strongly interacting normal 
state without any free parameters. Since we do not ex- 
pect the coefficients c± to depend strongly on tempera- 
ture, we can also use our method at low temperatures. 
Moreover, it is easily extended to other mass ratios and 
to the superfluid state ||. 

The above discussion shows that we effectively account 
for all fluctuations in the normal state. However, there is 
a second effect of particle-hole fluctuations which affects 
the superfluid state. Namely, the change in the coefficient 
a due to screening of the interspecies interaction [25j |. 
We account for this screening by replacing the two-body 
transition matrix T 2B with an effective transition matrix 
that includes the so-called bubble sum. We thus have 
that l/T off = 1/T 2B - 701(0,0), where 



MI(0,0) 



dk' N + (k') -iV_(k') 



(2tt) 3 2h' 



e-(k') 



(4) 



Since now a = -TiG^ (0,0) = -l/T ott + ?iS(0,0), the 
equation for the critical temperature, a = 0, becomes 
— 7iS(0, 0) = MI(0,0), where we use renormalized chem- 
ical potentials to include the selfenergy effects. If we ap- 
ply this procedure to the mass and population-balanced 
case we find T c = 0.18Tf and n(T c ) = 0.52ep, which is to 
be compared with the Monte Carlo results T c — 0.15Tp 
and (i(T c ) = 0.49eF [HI]. Moreover, we find for the 
mass-balanced tricritical point that /cbTc3 = 0.09eF+ and 
P C 3 = 0.25, which agrees with the experimental data Q. 
Returning to the 6 Li- 40 K mixture and using the same 
conditions for the tricritical point and the Lifshitz point 
as before, we find the result in Fig. lb. Note that the crit- 
ical temperatures are lowered by about a factor of three 
due to the screening and the selfenergy effects. Also the 
locations of the Lifshitz point and the tricritical point 



have drastically changed, since now we find T C 3 = 0.08Tp 
and P c3 = 0.47, while T L = 0.05T F and P L = -0.18. 

Conclusion. — We have considered the experimentally 
available 6 Li- 40 K mixture at unitarity, where we incorpo- 
rated selfenergy effects to reproduce at zero temperature 
the normal equation of state from Monte Carlo calcula- 
tions. By also including the effect of screening on the 
critical temperature, we have made quantitative predic- 
tions for the phase diagram, which contains both a Lif- 
shitz and a tricritical point. At weaker attractions the 
Lifshitz point remains present, although its temperature 
gets exponentially suppressed. Below it, various super- 
solid phases are competitive and a rich phase 
structure is expected. We hope that this Letter will stim- 
ulate new experimental research to supersolidity. 
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